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Abstract. Despite numerous efforts, the transition from Asymptotic Giant Branch (AGB) stars to Planetary 
Nebulae (PN) is a poorly understood phase of stellar evolution. We have therefore carried out interferometric 
(VLA) radio observations of a sample of hot post-AGB stars, selected on the basis of their optical and infrared 
properties. Ten sources out of the 16 observed were detected. This indicates that most of our targets are surrounded 
by a nebula where the ionization has already started. This definitively determines the evolutionary status of the 
selected sources and provides us with a unique sample of very young Planetary Nebulae (yPNe). A comparison 
with another sample of yPNe confirms our working hypothesis that our targets are indeed very young, probably 
just in the transition toward PN. Finding transition objects is extremely important as they can provide unique 
clues for a better understanding of this important phase of stellar evolution. 

Key words. Stars:AGB and post-AGB - Planetary Nebulae: general- radio continuum; stars. 



1. Introduction 

The formation and early evolution of Planetary Nebulae 
(PNe) is a quite obscure phase of stellar evolution. 
Recently HST has provided high-quality images of PNe 
which quite often reveal jets or bipolar structures whose 
origin it is difficult to explain. In particular, it is quite 
challenging to understand how the spherically symmet- 
ric Asymptotic Giant Branch (AGB) circumstellar shells 
transform themselves into the extremely asymmetric en- 
velopes observed in evolved PNe. 

Recent models propose different agents in shaping PNe 
such as high-speed, collimated outflows or jets that oper- 
ate during the early post-AGB evolutionary phase (Sahai 

& Trauger,[TnMI Lee & Sahai, EOOS- 
Very recently, high sensitivity, high resolution MERLIN 
observations of IRAS 18276-1431 in the OH maser lines 
1612 and 1667 MHz have made it possible, for the first 
time, to detect and measure the circumstellar magnetic 
field in a PPN (Bains et al., 2003). This result suggests 
that the magnetic field may play a fundamental role in 
the process that shapes the circumstellar envelopes in the 
post-AGB evolutionary phase. 

During the last few years many observational pro- 
grams have been devoted to recognizing new planetary 
and proto-planetary nebulae among unidentified IRAS 
sources with far infrared colors similar to those of known 
planetary nebulae (Pottasch et al.. 119881 van de Steene & 
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Pottasch. 119931 Garcia-Lario et al.. ll997a|l . The final aim 
of these studies is to fully understand the process of for- 
mation of Planetary Nebulae through the discovery and 
analysis of new transition objects in the short phase be- 
tween the end of the AGB and the onset of the ionization 
in the nebula. 

The number of known transition objects is extremely 
small, tens compared to more than a thousand PNe iden- 
tified in our Galaxy (Kwok. ri994|l . This is related to the 
short evolutionary timescales involved (a few thousand 
years) . 

As post AGB stars/proto-PNe evolve into early stages 
of Planetary Nebulae, part of the shell begins to be ion- 
ized by the central star. The post-AGB objects detected 
by IRAS extended from non variable OH-IR stars to M, K, 
G, F, A supergiant types. This sequence appears to rep- 
resent the evolution of the post-AGB stars towards hotter 
spectral types and we should also expect to find some B- 
type stars among the hottest post-AGB stars. There is a 
small sub-group of B-type stars, called BQ[ ] stars, defined 
as Be with forbidden emission lines, and Parthasarathy & 
Pottasch Ijl989|l . on the basis of their IR excess, recognized 
them as potential candidates to be new transition objects. 

BQ[ ] stars, however, are not a well defined group, and 
there is still a controversy on their evolutionary stage. 
Parthasarathy & Pottasch (1989) suggested that at least 
a fraction of them could be hot post-AGB stars in an 
early stage of their evolution as proto-PNe. The common 
characteristics of these kinds of star are: 



2 



Umana et al.: A search for very young Planetary Nebulae 



— numerous permitted and forbidden emission lines of 
several elements, in addition to the absorption lines 
typical of a Bl I-II star; 

— underabundance of carbon and metals; 

— high galactic latitude; 

— presence of PN-type detached cold circumstellar shells. 

These characteristics are typical for low-mass, highly 
evolved post-AGB stars. 

In some cases the central star may just have become hot 
enough to photoionize the circumstellar envelope ejected 
during the previous AGB phase, as the presence of low 
excitation nebular emission lines in the spectrum seems 
to indicate. 



2. Observations and Results 

2.1. Sample selection 

In order to find new very young Planetary Nebulae we 
have selected a sample of 16 hot post-AGB stars from the 
most recent compilations, namely Parthasarathy, QOO Sal: 
Conlon et al., and Parthasarathy et al., (jSODOa,). 

All the selected candidates have high galactic latitude, 
infrared excess, spectral type Bl I-II and nebular emission 
lines in their spectrum (Table In particular the last 
two requirements maximize the possibility of detecting a 
radio nebula. Moreover in one of the proposed targets, 
namely SAO 85766, the strong variations observed in the 
spectrum suggest that it has rapidly evolved in the last 50 
years (Parthasarathy et al.. 12000^)1 and its evolutionary 
state appears to be similar to that observed in the case of 
the hot post-AGB star SAO 244567, where a radio nebula 
has been detected (Umana et al.,[2004). 

The detection of radio nebulae associated to the se- 
lected targets would definitively establish the evolutionary 
status of this kind of objects, producing a unique sample 
of very young PNe. 

Successive higher resolution radio observations of the de- 
tected targets will allow us to study nebular morphologies 
in the very early phases of evolution of PNe and thus clar- 
ify at what stage of the AGB-PN evolution those asymme- 
tries are established , providing strong clues to the process 
shaping PNe. 

2.2. The VLA data 

The observations were carried out using the VLA ^ in two 
different runs, the first on June 8, 2001, from 06:06 UT 
to 11:33 UT, and the second on June 10, 2001 from 18:56 
UT to 20:24 UT. 

We observed our targets at 8.4 GHz (X-Band), using two 
independent 50 MHz bands. The choice of the observing 
frequency is a compromise between the sensitivity of the 
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instrument and the higher probability of observing in the 
optically thin region of the spectra. The array was in CnB 
configuration, providing a typical beam size of ~ 2.5". 

A typical observing cycle consisted of 10 min integra- 
tion time, preceded and followed by a 2-min observation of 
the phase calibrator. This basic sequence was repeated at 
least 3 times in order to improve the signal to noise ratio 
and thus to ensure the necessary sensitivity for the de- 
tection of weak sources. The flux density scale was deter- 
mined by observing 3C286 and 3C48, while phase calibra- 
tion was obtained from observations of phase calibrators 
close to our targets. 

The data processing was performed using the standard 
programs of the NRAO Astronomical Image Processing 
System (AIPS). To achieve the highest possible signal 
to noise ratio, the mapping process was performed using 
the natural weighting and the dirty map was CLEANed 
down as close as possible to the theoretical noise. To esti- 
mate the noise level in the maps we analyzed an area on 
the map, whose dimension is of the order of more than 100 
0gyjj , away from the phase center and free from evident ra- 
dio sources. Its consistency with the expected theoretical 
noise was also checked. 

Of 16 observed sources, we detected a total of ten. The 
source position and the flux density determination were 
obtained by fitting a gaussian brightness distribution 
(JMFIT). The uncertainty in the source position is a func- 
tion of the rms in the map, the resolution (^syn), and 
the peak flxrx density of the source in the image (i^peak), 
namely: 



rms 



2K 



peak 



"syn 

2^ 



(1) 



where is the signal to noise ratio for the source de- 
tection. For the worst case in our sample we derive a dpos 
of about 0'.'12. 

Our results are summarized in Table 1, where the radio 
flux density, with associated error (cr) and the rms of the 
map, are reported. The error associated to the flux density 
estimation is a combination of the rms in the map plus the 
amplitude calibration error. Therefore, a is derived from: 



O- = V' (rTOs)2 + ((Tcal^3.6cm)^ 



(2) 



where (Tcai is the error associated to the flux calibration, 
which is of the order of 2-3 %, and Fa.gcm is the measured 
flux density. At CnB resolution all the detected targets 

are point-like or slightly resolved radio sources, i.e. 

n ^ off 

^'source _ ^ • 

All the sources in our sample were checked for 20 cm 
emission using the NRAO VLA Sky Survey (NVSS, 
Condon et al., I1998|l . Only two of our detected targets 
were present in the survey, namely IRAS 18442 and 
IRAS 19336, with fluxes at 1.4 GHz of 8.4 ± 0.6 and 
8.8 ± 0.6 mJy respectively. 
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It is not easy to verify whether the non-detections are 
related to a larger distance or due to an intrinsic charac- 
teristic of the source such as the evolution of the ionization 
structure. In fact, distances to planetary nebulae are no- 
toriously very difficult to be estimate. 

Traditionally, distances are derived through statistical 
methods, which are based on a particular characteristic 
of the nebula, assumed to be invariant. For example, the 
Shklovsky method assumes that all PN have the same 
ionized masses (Cahn & Kaler [T!l7T|l . Since this is often 
not the case, statistical methods are considered not to be 
rehable CTerzian. lTO^ . 

Recently, Gauba et al. H2l)(J3fl modelled the near and 
far-infrared flux distribution of a small sample of hot 
post-AGB candidates, with IRAS colors similar to those 
of PNe. Their sample includes four of our targets, namely 
IRAS 17074-^1845, IRAS 17423-1755, IRAS 19157-0247 
and IRAS 19399+2312. Using the DUSTY code (Ivezic 
et al. 11999(1 . the authors modelled the circumstellar shells 
of their sample and derived various physical parameters, 
including an estimate of the distance. They derived a 
distance of 3.7, 3.7, 0.9 and 0.7 Kpc respectively to 
IRAS 17074+1845, IRAS 17423-1755 IRAS 19157-0247 
and IRAS 19399+2312, with an uncertainty of about 10% 
for each source, depending on the assumed core mass. 
Their results indicate, at least for a small fraction of our 
sample, that the lack of detectable radio flux should be 
related to intrinsic factor and not to the distance. As a 
matter of fact, the undetected sources IRAS 19157-0247 
and IRAS 19399+2312 have distances much smaller than 
the detected source IRAS 17423-1755. 

We stress again the difficulties in deriving distances 
to PNe, regardless of the used method, as it will lead to 
uncertainties by a typical factor of 2 or 3. As an exam- 
ple, the distance to IRAS 17423-1755 has been derived by 
different authors, using different methods, to range from 
8.3 kpc (Borkowski & Harrington. l2Q0l] l to 5.4 kpc (Riera 
et al. . 1200311 . Therefore, results presented by Gauba et al. 
H2UU3|I provide only an indication and should be taken 
with caution. 



3. Notes on individual sources 

To point out any characteristics common to the detected 
sources, besides selection criteria, we summarize, in 
the following, the physical properties of each source as 
derived from the literature. When available, a previous 
radio measurement is reported and compared with our 
result. All the targets in our sample are associated witn 
an IRAS source. When a discrepancy between radio, 
optical and IRAS coordinates is present, a comment is 
added at the end of each source section. 

IRAS 06556+1623 

When it was discovered, in 1932, this star showed 
some P Cygni line profiles in its spectrum. Ciatti et al. 



II1974|I pointed out the presence, in the spectrum, of sev- 
eral strong emission lines. The continuum emission was 
well represented by a BO. 5 star and the line ratios for 
[O i] , [O ii] , [O III] were not compatible with a unique set 
of density and temperature of the nebula. This suggests 
that the star is surrounded by an extended gaseous shell 
composed of regions of different temperatures, densities 
and velocities (Klutz & Swings. [l977|l . 
By comparing their spectrophotometric data with other 
data previously reported by other authors, Arkhipova & 
Ipatov H1982|l concluded that the volume and density of 
the gaseous shell had been increasing for the last 40 years. 
The more recent spectra are those obtained by Jaschek et 
al. ©96), which were collected between 1990 and 1993. 
Many lines are present, including several classic nebular 
lines such as [O ill] . No significant differences with previ- 
ous similar work are reported. Attempts to detect radio 
emission were made by Sistla & Seung f 1975|l at 8.4 GHz 
and by Altenhoff et al. H197()|l at 10.7 GHz, but both stud- 
ies the sensitivity was not high enoigh to detect the source, 
the rms being of the order of 5 mjy. 
IRAS 17381-1616 

The star was included in a sample of PN candidates se- 
lected from the IRAS Point Source Catalog on the basis of 
their far-IR colours fPreite-Martinez. [TM^ . IRAS 17381- 
1616 was included in a list of objects that were probably 
not PN but PN-related. Parthasarathy et al. ( 2000a) clas- 
sified it as a Bllbe star, with H/3 and H7 emission lines, 
suggesting that the object was a post-AGB star. 

The interferometer was pointed at LS 4331, which is 
believed to be the optical counterpart of IRAS 17381-1616. 
However, we did not find the radio source at the expected 
optical coordinates, but we detected a 1.42 ± 0.05 mJy 
radio source at a radio position of a = 17 : 41 : 00.05 and 

5 — —16 : 18 : 12.45, corresponding to a shift of 2.16" in 
right ascension and about 16" in declination, or about 6.5 

^syn ■ 

However the radio positions are shifted by only Aa — 
2" 16 and AS = 1'.'15 with respect to the IRAS coordinates, 
well inside the IRAS ellipse error (34" x 8" ) . We can thus 
conclude that our radio source is associated to the IRAS 
source. 

IRAS 17423-1755 

The source was believed to be a massive population I 
B-type star, but Bobrowsky et al. ((1995y) recognized it as 
a new planetary nebula. 

High resolution spectroscopy and optical images re- 
veal a complex bipolar morphology: a dense central region 
(disk) from which large-scale ffows collimate into bipolar 
jets (Bobrowsky et al. 119951 Riera et al.. 119951 Borkowski 

6 Harrington, 200T|). 

Parthasarathy et al. I|2()()()a|l found Balmer lines having 
a P Cygni profile, blue-shifted by — 400kms~^ and from 
the analysis of H^ profiles Sanchez Contreras & Sahai 
1)2001(1 pointed out the presence of two different winds, 
one of these probably associated to the post-AGB outflow 
prior to its interaction with the AGB wind. 
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Table 1. Programme Stars 



Name 


IRAS 


QJ2000 
1 




<5j2000 

1 


Q^radio 




f^radio 


V. 


S.T. 


-fl^ .6cm 

mJy 


mJy 


rms 
mJy 


OY Gem 


06556-1-1623 


06 


58 


30 


3 


+16 19 25 


06 


58 


30 


3 


+16 19 26.2 


t -1 ffQl 

11.5^^ 




0.55 


0.03 


0.02 


Hen 3-1347 


17074-1845 


17 


10 


24 


1 


-18 49 01 












9.5'^ 








0.03 




17203-1534 


17 


23 


11 


8 


-15 37 15 












12.5^'*^ 








0.03 


LS 4331 


17381-1616 


17 


40 


59 


9 


-16 17 58 


17 


41 


00 


1 


-16 18 12 


12.2*''^ 


Bile'-*' 


1.42 


0.05 


0.03 


Hen 3-1475 


17423-1755 


17 


45 


14 


2 


-17 56 47 


17 


45 


14 


1 


-17 56 45 


11.0^ ' 


Be'-''' 


0.26 


0.03 


0.03 


SAO 209306 


17460-3114 


17 


49 


16 


5 


-31 15 18 


17 


49 


16 


5 


-31 15 18 


7.9^ 


BOIIe^ ' 


1.29 


0.05 


0.04 


V886 Her 


18062-f2410 


1 Q 

if5 




on 


1 


+24 10 43 


1 Q 

io 




on 


1 


+24 10 43 


11.6*''* 


Blllle'- ' 


1.46 


0.05 


0.03 


LS 63 


18371-3159 


18 


40 


21 


6 


-31 56 49 


18 


40 


22 





-31 56 49 


ILg**** 


T->-1 T (A) 

Bile' ' 


0.62 


0.03 


0.03 


T,S 51 1 9 


1 837Q-1 707 


18 


40 


48 


5 


-17 04 36 




















0.03 




18435-0052 


18 


46 


06 


9 


-00 48 55 












ii.ow 


B2//(*' 






0.03 


BD -11 4747 


18442-1144 


18 


47 


04 


8 


-11 41 02 


18 


47 


04 





-11 41 12 


9.3(2) 




19.21 


0.60 


0.05 


LS IV -02 29 


19157-0247 


19 


18 


22 


5 


-02 42 09 












II.5W 


Bim'*' 






0.03 




19336-0400 


19 


36 


17 


5 


-03 53 24 


19 


36 


17 


5 


-03 53 25 


12.5'''* 


Bl/e(^) 


9.74 


0.29 


0.03 


LS H -1-23 17 


19399+2312 


19 


42 


05 


5 


+23 18 59 












IO.4W 


Bl/J/W 






0.03 


LS IV -12 111 


19590-1249 


20 


01 


49 


8 


-12 41 17 


20 


01 


49 


8 


-12 41 17 


11.3(''> 


Bl/e(4) 


2.76 


0.08 


0.03 


LS II +34 26 


20462+3416 


20 


48 


16 


6 


+34 27 24 


20 


48 


16 


6 


+34 27 24 


ll.l(^) 




0.42 


0.03 


0.03 



1) Optical coordinates are FK5 2000.0 as reported by SIMBAD database. 

2) Parthasarathv. I1993al 

3) Parthasarathv. I1993bl 

4) Parthasarathy et al.. l2()()()al 



Knapp et al. H1995(l detected it at 8.4 GHz with the 
Very Large Array (NRAO) and found the flux to be 
0.30 ±0.04 mJy, with an excitation parameter of 1.8, com- 
patible with a B3 spectral type. Our flux determination is 
in agreement, within the errors, with previous results. 
IRAS 17460-3114 

Crampton (,197Tjl classified it as an 07.5 star in a 
Hii region; this was confirmed by Dachs et al. fl9S2|. 
Gies and Mason et al. ijTO^ found a 08V spec- 

tral type. On the basis of its far-IR colors similar to 
PNe and its spectral type Parthasarathy H1993afl classi- 
fied IRAS 17460-3114 as a hot post-AGB star. 

Radio measurements were made by Ratag et al. I)1990f) . 
who found a flux of 0.8 mJy at 4.8 GHz using the VLA. 
We detect this source with a flux density at 8.4 GHz of 
1.29 ± 0.05 mJy. Combining our results with data from 
Ratag et al. (1990), we derive a spectral index a — 0.85, 
consistent with an optically thick nebula, at least between 
4.8 and 8.4 GHz as we expected from a very young, and 
thus dense, PN. 
IRAS 18062+2410 

According to the HDE catalogue its spectral type was 
A5 in 1940, but Downes & Keyes (,1988.) classified it as a 
Be star. Arkhipova et al. H1999|l and Partasarathy et al. 
Ij2000b|l analyzed the high resolution spectrum of this star 
and found several permitted and forbidden emission lines: 
the presence of [N 11] and [S 11] lines indicate the presence 
of a low excitation nebula. Rapid photometric variability 
has also been detected (Arkhipova et al.. I2001a|l which. 



together the observed spectral variations, appears to in- 
dicate that the star has been evolving rapidly in the last 
150 years. 

IRAS 18371-3159 

Preite-Martinez 1)1988(1 selected this star from the 
IRAS PSC as a planetary nebula candidate on the basis 
of its far-IR colors. Parthasarathy et al. H2000a|l found its 
spectral type to be Bllabe, which combined with its high 
galactic latitude, far-IR colors similar to PN, and Balmer 
line emission, makes it a post-AGB candidate. 
IRAS 18442-1144 

This source is very near to BD-11 4747 and possi- 
bly constitutes a binary system or is identifiable with it. 
Corporon & Lagrange ( 19^^ did not find any evidence 
of binariety. Parthasarathy et al. ( 2000a ) found an A3V 
spectral type for BD-11 4747, but they could not observe 
the optical counterpart of IRAS 18442-1144 because it was 
probably too faint to detect. On the basis of the presence 
of [Nil] fines in the spectrum IRAS 18442-1144 may be a 
post-AGB star or a low excitation PN. 

We found no source at the coordinates of BD- 
11 4747, which is reported as the optical counterpart 
of IRAS 18442-1144. A distinct object was detected at 
i?A = 18 : 47 : 04 and DEC = -11 : 41 : 12, with a flux 
density of 19.21 ± 0.06 mJy. 

The radio positions are shifted by only Aa — 3" and 
A6 = 1" with respect to the IRAS coordinates, well inside 
the IRAS ellipse error (30" x 7"). We can thus conclude 
that the radio source is associated to IRAS 18442-1144. 
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If our data are combined with those from NVSS it appears 
that the nebula around IRAS 18442-1144 is opticaUy thick 
between 1.4 and 8.4 GHz. 
IRAS 19336-0400 

Downes & Keyes 1)1988(1 found [N ii] and [S ii] hues in 
the spectrum of IRAS 19336-0400, suggesting the presence 
of a nebula. Parthasarathy et al. (^2000aj) classified it as a 
Bl supergiant with Balmer lines in emission, but without 
[O III] line, thus suggesting it is a very young and low ex- 
citation planetary nebula. IRAS 19336-0400 was detected 
by van de Steene & Pottash ifTO^ with the WRST at 
4.8 GHz as a 10. 5± 1.1 mJy radio source. The source is 
also present in the NVSS with a flux density at 1.4 GHz of 
8.8 ± 0.6 mJy. We measured a 8.4 GHz radio flux density 
of 9.74 ± 0.3 mJy 

Combining our data with those from the literature we de- 
rive a quite flat radio spectrum, indicating a nebula which 
is optically thin between 1.4 and 8.4 GHz. 
IRAS 19590-1249 

Ultraviolet, optical and infrared observations of 
IRAS 19590-1249 (McCausland et al. lTM^ Conlon et al. 
119931 suggest that this BO supergiant is in the post-AGB 
evolutionary stage, evolving into a low excitation plane- 
tary nebula. By comparing the atmospheric parameters of 
the stars with theoretical post-AGB evolutionary tracks 
a central mass of 0.67 Mq has been estimated (Conlon et 
al. I1993|l . implying that further evolution can occur on a 
timescale not more than hundreds of years. Photometric 
variability has been reported (Arkhipova et al., 12002(1 . 
with a pattern similar to what has been observed in other 
two of our detected targets, namely IRAS 18062-f2410 
and IRAS 204624-3416, which appears to be attributable 
to their common evolutionary status. 

IRAS 20462-1-3416 

Turner & Drilling 1(1984(1 classifled this star as a B1.5 
supergiant, but, on the basis of its high galactic latitude, 
far-IR colours similar to PN, dust shell parameters and 
flux distribution, Parthasarathy | jl993b() . concluded that 
it was a low mass post-AGB star, with a detected cold 
dust sheU. The optical spectrum of IRAS 20462-f3416 
consists of the spectrum of a B1.5 supergiant plus the 
emission spectrum of a low-excitation nebula, with no 
[Oiii] lines (Parthasaraty, [1993b[ Smith & Lambert, 
I1994(l . Strong spectral variations were observed (Smith 
and Lambert, 119941 Garcia-Lario et al., 1997b^ . which 
are interpreted in terms of mass-loss episodes. Evidence 
of mass-loss has been reported also by Arrieta & Torres- 
Peimbert (.,20021 from the analysis of P Cyg profiles 
observed in both UV and optical spectra. As in other hot 
post-AGB star, the star exhibited rapid, irregular light 
variation, that can be attributed to a variable stellar wind 
and/or to some kind of stellar pulsations (Arkhipova et 
aL. 12001 bll . 

We want to stress that only few of our targets have 
been studied in detail. In particular, multi-epoch pho- 
tometric and spectroscopic studies are available only for 



IRAS 06556-H1623, IRAS 17423-1755, IRAS 18062-^2410, 
IRAS 19590-1249 and IRAS 20462-^3416. Such studies in- 
dicated a strong variability in both continum and spectral 
lines. 

We may conclude that, besides selection criteria and the 
observed free-free emission, a possible common character- 
istic to our sources is strong spectral variability which may 
indicate fast evolution through episodic mass-loss events. 

4. Nebular Characteristics 

Since the aim of this work is to find new objects in the 
very initial stage of PNe evolution, we need to check if 
our detected objects are indeed young PNe. We will first 
determine the radio characteristics of our sample and then 
we will check how they compare with those of young PNe 
as derived from previous, more extended radio surveys. 

4.1. Brightness Temperature 

/^From measurements of the radio continuum of PNe it is 
possible to determine a distance-independent parameter, 
namely the brightness temperature (Tb), which is defined 



(f F 
Tr = (— ) — 



(3) 



where F^, is the flux density at frequency v, and 9 is the 
angular diameter of the nebula. For ly = 8.4 GHz Eq.JJl 
becomes: 



Tr 



25(- 



-)■ 



-2 -^8.4 GHz 

mJy 



(4) 



At the resolution of our observations 4 detected targets are 
point-like. In these cases we can derive only a lower limit 
for the brightness temperature, assuming for the sources 
an angular diameter < 2". For the rest of the target the 
angular size has been derived by a Gaussian fitting of 
the source in the map (see Sect 2.2) and by performing 
a Gaussian deconvolution of the fitted source size. 

If we indicate with 6'syn the FWHM of the point-source 
response of the interferometer and with the fitted 
FWHM size, the resultant deconvolved FWHM size {9) 
is: 



^J{9Gr {9syny 



(5) 



which, for resolved sources, can be significatively smaller 
than the nominal resolution of the interferometer. 

The obtained angular sizes were then corrected for the 
appropriate factor (7) following van Hoof 1(2000(1 . 



7(/3) 



0.36 



1 + 0.7994/32 



1.1106 



(6) 



where (3 represents the ratio between 9g and 6'syn. 
The results are summarized in Table HJ As the nebula 
expands, as a consequence of PN evolution, its surface 
brightness is expected to decrease. Thus Tb is usually con- 
sidered an age indicator, as high Tr should characterize 
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Table 2. Summary of nebular characteristics of detected targets. Angular diameters are corrected following van Hoof 
H2000(l : The mean emission measure (EM) has been calculated from the formula of Terzian & Dickey H1973|l : IRE are 
derived following Pottasch H1984fl 



IRAS ID 


Diameter 


Tb 


EM 


IRE 


-^dust 




[arcsec] 


[K] 


[lO^'cm^'^pc] 




[K] 


06556+1623 


2.1 


2.3 


6.3 


194 


181 


17381-1616 


< 2.0 


> 8.9 


> 18.8 


31 


218 


17423-1755 


< 2.0 


> 1.6 


> 3.4 


2984 


130 


17460-3114 


1.1 


27 


56.6 


248 


204 


18062-h2410 


< 2.0 


> 9.1 


> 19.3 


106 


253 


18371-3159 


< 2.0 


> 3.9 


> 8.2 


187 


187 


18442-1144 


1.8 


148 


314.8 


11 


129 


19336-0400 


1.5 


108 


229.8 


14 


187 


19590-1249 


1.9 


19.1 


40.6 


21 


134 


20462-h3416 


2.2 


2.2 


4.6 


186 


114 



compact, very young PNe. Tb measured for large samples 
of PNe are in the range between 10000 and 0.1 K (Kwok, 
I199U|I . Our results indeed indicate very low Tb values for 
our sample. This, however, can be still in agreement with 
the hypothesis of young age if one assumes that at the ob- 
serving frequency the optical depth of the nebula is very 
small. 

Among the targets of our sample there are three 
sources, namely IRAS 17460-3114, IRAS 18442-1144 and 
IRAS 19336-0400, whose radio spectrum can be inferred 
on the basis of data from the literature and the results 
presented in this paper (see Sect 3) . It can be argued that 
the first two sources appear to be optically thick, at least 
between 20 and 3.6 cm. Unfortunately, the available spec- 
tral information is too poor and a two-point spectrum, 
in one case obtained very far apart in frequency, is not 
sufficient to draw any conclusion on the spectral index. 

4.2. Emission l\/leasure 

Another physical parameter whose high value should char- 
acterize a young PN is the electron density which is ex- 
pected to decrease as the nebula ages and expands. 
The mean emission measure (< EM >), as derived from 
radio observations, is an averaged property of the ionized 
nebula and, therefore, may better represent the overall 
density structure of PNe. 

If the source is optically thin we can use the observed 
radio flux density plus its angular size to derive the mean 
emission measure as follows (Terzian & Dickey, I1973|l : 



< EM >= 



EMdn 5.3 X 10^i^8.4 GHz 



n 



612 



(7) 



where i^g 4 qhz is the measured radio flux density, in mJy, 
at 8.4 GHz, and 9, in arcsec, the angular dimension of 
the radio emitting region at that frequency. The mean 
emission measure is in cm~^pc. 

Young PNe should have emission measures of the order of 
10*^ - lO^cm-^pc (Terzian & Dickev. .1973. Kwok et al.. 



11981(1 and, as reported in Tabled our sample, in general, 
appears to have a more diluted nebula. 

4.3. Infrared Excess 

A young planetary nebula is still surrounded by the rem- 
nant of the dusty envelope of its AGB progenitor and 
emits most of its energy in the far infrared as thermal radi- 
ation from heated dust grains. The infrared excess (IRE) 
is defined as the ratio of the observed total far infrared 
flux (-Fir) over the expected total infrared flux. On the 
hypothesis that the far infrared flux is due to thermal 
emission from dusty grains heated by Lya photons, this 
ratio is unity. In young and compact planetary nebulae 
dust heating by direct absorption of starlight is impor- 
tant, and therefore this ratio can be much larger than 
unity. For this reason the infrared excess is often used as 
an age indicator for PNe. Since the expected total infrared 
flux can be expressed in terms of optically thin radio flux 
density, we may derive IRE from the formula given by 
Pottasch H1984|l . In particular, as the nebulae associated 
to our targets are compact, the infrared excess is derived 
in the high density approximation, i.e.: 



IRE = 1.07- 



Fir 



4 GHz 



(8) 



where Fir, in 10^^'' Wm^^, is the total infrared flux. 
For consistent comparison with other similar studies, 
has been derived by fitting a Planck curve over the IRAS 
bands, namely from 12 to 100 fim (Aaquist & Kwok, 
11990(1 . Our sample is characterized by IRE values much 
higher than unity. As a by-product of this method we also 
derive the dust color temperature. A summary of the neb- 
ular characteristics of our sample is reported in Table [5] 

5. Comparison with another sample 

In order to check if our sample consists of young PNe, 
we will compare the radio and infrared properties of the 
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Fig. 1. Distribution of brightness temperature (Tb) for 
the AK91 sample (dashed hne) and the sample of the 
present paper (full line). Tb for AK91 were calculated from 
Eq.Q and using radio measurements obtained at 15 GHz. 
Four of our targets have only lower limits for Tb 

targets belonging to our sample with results obtained in 
similar studies. In particular, in the following, we will con- 
sider the work of Aaquist & Kwok (,1991t AK91) carried 
out with the VLA at 15 GHz on a sample of young PNe 
selected on the basis of their compact radio morphology. 
The targets of AK91 were also observed at 5 GHz, but 
we prefer to compare our results with those obtained at 
15 GHz since in both cases optical depth effects should 
not be important. 

All the selected targets in AK91 have high brightness tem- 
peratures, infrared excess (IRE) much higher than unity 
and dust temperature higher than the typical value ob- 
served in more evolved nebulae, which is of the order of 
100 K (Pottasch et al., 1984). All these properties are 
consistent with the hypothesis that the sample consists of 
very young PNe. 

In order to compare the physical properties of the 
nebulae belonging to AK91 with our sample we plot in 
Figs, n [3 and 01 the brightness temperature (Tb), the 
emission measure (EM) and the infrared excess (IRE) of 
both samples. Those quantities were re-calculated from 
the published radio measurements using the same formu- 
las as for our sample. 

It is evident that for the AK91 sample Tb and EM 
are systematically higher than for our sample, and this 
seems to indicate that our sample indeed consists of more 
evolved PNe. 

On the contrary, the infrared excess of our sample, which 
has values systematically higher than those reported by 
AK91, indicates a particularly young sample of PNe. This 
apparent contradiction is further complicated by the fact 
that the infrared properties of both samples are quite sim- 
ilar, as is evident from an inspection of the IRAS color- 
color diagram (Fig. 2]), where sources belonging to dif- 
ferent samples occupy the same region of the diagram. 
This region is also shared with SA0244567, the youngest 



Fig. 2. Distribution of mean emission measure (EM) of 
the AK91 sample (dashed line) and of the present sample 
(full line). Four of our targets have only lower limits for 
EM 
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Fig. 3. The radio flux density versus the far infrared flux 
for the present sample (empty triangle) and for the AK91 
sample (plus signs) 

known PN, whose evolution appears to be quite rapid, 
since it has become ionized only within the past 20 years 
(Parthasarathy et al.. I1993|l . 

In the diagram the infrared colors are defined as usual 
(van der Veen & Habing. 1 1988)1 : 

[12] -[25] = -2.5 log |H (9) 
where F12 and F25 are the fluxes measured at 12 and 25 

Moreover, dust temperatures of both samples have quite 
similar distributions (Fig.|SJ), strengthening the conclusion 
that both samples have similar dust characteristics. 
Concerning the possible causes of lower Tb and EM of our 
sample we may consider whether this could be related to 
a systematic effect due to the different spatial resolution 
used in the two surveys, as both Tb and EM depend on 
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the angular source size (oc 6'"^). 

A factor 20 in spatial resolution yields a two orders of 
magnitude difference in both Tb and EM, forcing the two 
distributions to be more similar than what is shown in 
Figs, n] and |21 However, only 4 out of 10 detected sources 
were not resolved and thus without a good angular size 
estimate. Therefore, while this effect cannot be ruled out, 
at least for these 4 sources, we will consider other possi- 
bilities which may explain our results. 

Volk & Kwok H1985|l studied the dynamical evolu- 
tion of a nebula in the framework of the interacting wind 
model. They considered all the possible factors that can 
affect the dynamics of the ejecta, including radiative cool- 
ing, variation of mass-loss of the stellar wind as well as 
evolution of the central object. 

The authors also reproduced, in the context of their dy- 
namical model, the time evolution of the radio flux density, 
assuming that all the radio emission arises only from the 
ionized part of the shell. 

The evolution of the radio flux density with the age of 
the nebula is outlined in Fig. adapted from Fig. 15 of 
Volk & Kwok 1)1985(1 . where three distinct phases of radio 
emission have been indicated. 

In particular, during the first phase, when the radio emit- 
ting region is bounded by the ionized front, the radio flux 
increases with the radius of the ionized nebula, which is 
growing with the age of the nebula. This is strictly related 
to the rate of Lyman continuum photons which increases 
with the central star evolution (Kwok, 1985). 

In this context we may try to solve the apparent in- 
consistency between the typical values of Tb and EM, as 
derived for our sample, which indicate a more evolved sam- 
ple, and the IRE values, which, on the contrary, appear 
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Fig. 5. Distribution of dust temperature of the AK91 sam- 
ple (dashed line) and of the present sample (full line) 

to be consistent with a population of very young PNe. 
If we consider the functional dependence of Tb, EM and 
IRE on the radio flux density (Fradio), the first two being 
proportional to Fradio and IRE to -Fj-adio' '^^'^ contradictory 
result can be explained if we assume that, on average, our 
sample is characterized by a very low radio flux compared 
to what was measured in the AK91 sample. 
We thus conclude that the present sample is, indeed, very 
young and that the ionization in the circumstellar nebula 
has just started. 

Therefore, only a small fraction of the nebula is ionized 
and in the evolutionary scheme, reproduced in Fig.|Sl they 
occupy the very initial portion of the ionization bounded 
phase, corresponding to a very low value of radio flux den- 
sity. On the contrary, the AK91 sample, even if consisting 
of young PNe, occupies a more advanced portion of the 
same evolutionary phase, corresponding to a much higher 
radio flux. 

To further test our hypothesis, we calculate the ionized 
mass of the nebula from the observed radio flux, following 
Pottasch nMl|): 

Mi = 5.11 X 10~''^(F8.4 61^5)1/2 (;lo) 

where it is assumed that the nebula has a filling factor 
e ~ 0.6 and that 10% of its mass consists of helium atoms. 
In the formula F^,^ is in mJy, 9 is in arcsec, d in kpc and 
Mi in Mq- If we assume a standard distance of 1 kpc for 
the sources detected in our survey, we obtain values for the 

ionized mass of the nebula in the range 6 x 10"^ < -M^ < 

- Mq - 

5.4 X 10~^, which are much lower that the typical value 
of 0.2Mq obtained for more evolved nebulae (Pottasch, 
I1984f) . If we assume that all the nebulae are at a standard 
distance of 10 kpc, we obtain values 2 x 10^^ < -M^ < 

- Mq - 

17 x 10~^. These are to be considered upper limits. 

A radio nebula was recently detected around 
SA0244567 (Umana et al.. l2004|l and preliminary results 
were presented by Trigilio et al. (lUUT). A Tb = 700 K 
and a EM = 1.5 x lO'^cm^'^pc place SA0244567 between 
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Adapted from Volk & Kwok H1985|l . The curve is obtained 
assuming a nebula at a distance of 1 Kpc, with a core- 
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our sample and AK91 and this implies, in the framework 
of our hypothesis, that this young PN possesses a radio 
nebula slightly more evolved than those in our sample. 



quite obscure evolutionary phase. 

To further test our results, we made a comparison with 
another sample of PNe assumed to be young. We preferred 
the sample studied by Aaquist & Kwok H1991|l . as both 
samples were observed at high frequencies and any pos- 
sible opacity effect can be neglected. By comparing the 
nebular characteristics of both samples, as derived from 
radio measurements, we conclude that our sample com- 
prises extremely young PNe. The apparent inconsistency 
between radio and infrared properties can be explained 
by assuming that in our nebulae the ionization has just 
started, while in the nebulae from AK91 the ionization 
structure is more evolved. This translates into a difference 
in the radio emission, which is, on average, less intense for 
our sample. This result is strengthened by the fact that 
we derive values of ionized mass of the nebula which are 
much lower than those usually derived for more evolved 
nebulae. 

Successive multi-frequency and high spatial resolution ra- 
dio observations will allow us to fully characterize the ra- 
dio properties of these new objects, providing important 
clues to better understand the formation and shaping of 
PN in the early stages of their evolution. 
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SIMBAD database, operated at CDS, Strasbourg, France. 



6. Summary 

The formation and early evolution of Planetary Nebulae is 
an unclear phase of stellar evolution. It is still not well un- 
derstood how the circumstcUar envelopes, which form dur- 
ing the AGB phase, transform themselves into the com- 
plex morphologies often observed in evolved PNe (Sahai, 
Kwok et al., '20025. In order to understand the ori- 
gin of such morphologies, it is very important to clarify at 
what stage of PN evolution the nebula becomes asymmet- 
ric or bipolar. Moreover, recently some objects have been 
identified whose evolution towards the PN phase occurs on 
timescales twenty times shorter than what is foreseen by 
current evolutionary models fParthasarathv. l^OOO|l . In this 
context it appears to be quite important to find objects 
in transition towards the PN phase, which can provide 
strong constraints on determining which process shapes 
the nebula and how this influences the nebular evolution. 

In order to improve our knowledge in this field we 
have observed a sample of 16 hot post- AGB stars, selected 
on the basis of their optical and infrared properties. A 
common feature of our targets is an infrared excess, sig- 
nature of mass-loss that occurred during the AGB phase. 
The main result of the present work is the detection of 
free-free radio emission in ten of the observed sources. 
This indicates that ionization has already started in their 
circumstellar material and that the sources can be con- 
sidered to be in a very early stage of Planetary Nebulae 
evolution. Therefore, the detected sources constitute a 
unique sample to be studied to for shedding light on this 
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